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ABSTRACT 

We present HST /WFC3 narrowband imaging of the Ha emission in a sample of eight 
gravitationally-lensed galaxies at z = 1 — 1.5. The magnification caused by the fore- 
ground clusters enables us to obtain a median source plane spatial resolution of 360pc, 
as well as providing magnifications in flux ranging from ~ 10 x to ~ 50 x . This enables 
us to identify resolved star-forming Hll regions at this epoch and therefore study their 
Ha luminosity distributions for comparisons with equivalent samples at z ~ 2 and in 
the local Universe. We find evolution in the both luminosity and surface brightness of 
Hn regions with redshift. The distribution of clump properties can be quantified with 
an Hn region luminosity function, which can be fit by a power law with an exponential 
break at some cut-off, and we find that the cut-off evolves with redshift. We therefore 
conclude that 'clumpy' galaxies are seen at high redshift because of the evolution of 
the cut-off mass; the galaxies themselves follow similar scaling relations to those at 
z = 0, but their Hn regions are larger and brighter and thus appear as clumps which 
dominate the morphology of the galaxy. A simple theoretical argument based on gas 
collapsing on scales of the Jeans mass in a marginally unstable disk shows that the 
clumpy morphologies of high-z galaxies are driven by the competing effects of higher 
gas fractions causing perturbations on larger scales, partially compensated by higher 
epicyclic frequencies which stabilise the disk. 
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1 INTRODUCTION 

Observations of star-forming galaxies at high- 2 have shown 
that a significant fraction of the population have tur- 
bulent, clumpy, rotating disks with clump masses of 
~ 1O 8 " 9 M , a factor of ~ lOOx the typical Gi- 



ant Molecular Cloud (CMC) locally (e.g. Cowie et al.l 
19951; lElmegreen et al.ll200l; l Elmegreen fc Elmegreerj|2005l ; 



Forster Schreiber et all 120091 ; lElmegreen et all 120091 '). The 
clumps are thou ght to form from gravitational instabilities 
in gas-rich disks (lElmegreen et aL I I20071 . 120091 ; iGenzel et al.l 
120081 ; iBournaud et alj|20ich . 

Some recent numerical simulations have suggested that 
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the majority of massive, high-z galaxies accrete their gas 
via 'cold flows,' in which the gas is accreted smoothly along 
filaments. These cold flows are less disruptive than a ma- 
jor merger, and hence offer a route to maintain marginally 
stable disks (Toomre parameter Q ~ 1) without disrupt- 
ing the structure and dynamics. Cold-flow accretion is ex- 
pected to be a dominate mode of mass assembly above z ~ 1, 
and thus acc ounts for the ubiquity of large clumps at high 
redshift (e.g. Bournaud fc ElmegreenlboOS ; IBournaud et al.l 
l201ll ; lDekel et al.ll2009l l~ 

In this picture, the clumps are considered to be tran- 
sient features, forming in marginally unstable disks at high- 
z and fed by smooth accretion of gas onto the galaxy. 
Clumpy galaxies therefore represent a phase in the evolu- 
tion of present-day spiral disks. 
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Figure 1. HST/ACS and WFC3 three-colour images of the observed clusters with the critical line at the redshift of the target arc 
overlaid, showing the positions of the target arcs. The arcs are contained within the white dashed boxes which denote the regions 
extracted in Figure [2] 
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There is a need to test the internal physical properties of 
the interstellar medium (ISM) observationally, to determine 
whether the clumps are scaled-up analogues of local Hll re- 
gions or represent a different 'mode' of star formation, and 
whether they can explain the strong evolution of star for- 
mation rate density with redshift. However, sufficient spatial 
resolution is required to resolve the ISM on the scales of star- 
forming regions. Even with the use of adaptive optics, spa- 
tially resolved studies of high-redshift galaxies to date have 
been limited to a resolution of ~ 1 .5kpc (e.g. iGenzel et al.l 
l200d : iForster Schreiber et alj|200sh : using the Hubble Space 
Telescope (HST), only the large st starburst complexes can 
be resolved, on scales of ~ lkpc (|Elmegreen et al.ll2007T l. On 
these scales, it is p ossible to probe the dynamics of galaxies 
on large scales, and lGenzel et all (|201ll ) found evidence that 
Q < 1 in the regions of galaxies where clumps are found, 
lending observational support to the theory that the clumps 
form from internal gravitational instabilities. In order to 
study the clumps in detail, we need to resolve high-redshift 
disks on the scales of individual star-forming regions; in the 
local universe, this is ~ lOOpc. 

The required spatial resolution can currently only be 
achieved by exploiting gravitational lensing. By targeting 
galaxies that lie behind foreground cluster lenses, it is pos- 
sible to benefit from linear magni fication factors (along 
one direction) of up t o 50x (e.g. ISwinbank et all 120071 
120091 ; IJones et alj l201Ch . and to isolate Hn regions of or- 
der ~ lOOpc out to z ~ 5 |Swinbank et alj |2009| ). Re- 
gions were found with star formation surfac e densities Ssfr 
~ 100 x higher than t hose found locally (|Swinbank et al.l 
120091 ; |j ones et alJl201Ch . These regions of dense star forma- 
tion are comparable to the most inten sely star-forming in - 
teracting systems in the local Universe (jBastian et al.|[2006f ) , 
yet appear to be ubiquitous in non-interacting galaxies at 
high redshift. 

It is not known what driv es these regions of i ntense star 
formation at high-z, although I Jones et al.l (|2010h suggest a 
combination of higher gas density, increased star formation 
efficiency and shorter star-formation timescales. In addition, 
their data give the appearance of a bimodal distribution of 
Hll region surface brightnesses, although there is no known 
physical process that might drive this. In order to under- 
stand this result further, we require a sample at interme- 
diate redshift (2 ~ 1 — 1.5) with which we can probe the 
evolution of star formation density with redshift at higher 
sensitivity so that regions comparable to those at z = are 
detectable. 

Previous work on high-z clumps has made use of 
Integral Field Units such as Keck /OSI RIS (| Jones et al 



2010[ : IWisnioski et alj2012h , Gemini /NIFS dSwinbank et al. 
20091 ) and VLT/SINFONI (|Forster Schreiber et all l2009h . 



These allow detailed mapping of the nebular emission lines, 
but at lower sensitivity than is achievable with imaging. 
An alternative means of identifying star-forming regions 
with high sensitivity is to take imaging through narrow- 
band filters. The Wide Field Camera 3 (WFC3) on the 
HST presents an opportunity to study the star formation in 
galaxies at z ~ 1 and z ~ 1.5, as there are narrowband fil- 
ters available which correspond to the wavelength of the Ho 
emission line at these redshifts. Combining the sensitivity 
and high spatial resolution of HST /WFC3 with the mag- 
nification afforded by gravitational lensing by foreground 



clusters, we can map the internal star formation distribu- 
tion and so identify the frequency and properties of giant 
Hll regions. 

In this paper, we therefore study the star formation 
morphologies of eight galaxies at z ~ 1 — 1.5. We present 
the sample in fJ3] present the properties of the galaxies and 
their star-forming clumps in ij3] discuss the implications in 
i|4]and present our conclusions in fj5] Throughout, we adopt 
a ACDM cosmology with H = 70kms~ 1 Mpe -1 , fi A = 0.7 
and O m = 0.3. Star-formation rates are calculated f r om Ha 
luminosity L u_ a using the pre scription of iKennicuttl (|l998aT ) 
adjusted to a lChabrierj <|2003h IMF. 



2 SAMPLE AND OBSERVATIONS 

Our sample comprises eight lensed galaxies, each with 
spectroscopically-confirmed redshifts between 1 < z < 
1.5 such that the Ha emission line falls within the high- 
transmission region of the narrowband filters on WFC3. The 
associated cluster lenses are massive s ystems from the X- 
ray s e lected BCS and MACS samples (|Ebeling et"allll998l . 
120011 . 120071 . 12010? ) with well-constrained mass models (see 
references in Table [TJ , so that the effects of lensing can be 
accounted for. 

The positions and properties of the sample are given in 
Table [T] We observed each target in the narrowband filter 
covering Ha for a typical exposure time of 6ks (2 orbits), 
using a 3- or 4-point linear dithering pattern of ±5 arcsecs 
in both directions to improve the detection and removal of 
cosmic ray hits and bad pixels. At the same time, three of 
the targets (MACS J0947, MACS J0159 and MACS J1133), 
which did not have WFC3 data in the archive, were ob- 
served in the corresponding broadband filter using the same 
sequence of observations as their corresponding narrowband 
data, for a total of 3ks (1 orbit). The narrowband data 
and new broadband observations were obtained in Cycle 18 
under Program 12197 (PLRichard), with the exception of 
Abell 2390, for which the broadband and narrowband data 
were taken in Cycle 17 under Program 11678 (PLRigby). 
The remaining broadband data were obtained under Cycle 
17 Program 11591 (PLKneib) or Cycle 18 Program 12065-9 
(PLPostman) as indicated in the notes to Table [T] 

All of t he WFC3 data were redu ced using the multidriz- 
zle software (|Koekemoer et al"1l2002h under PyRAF to perform 
a cosmic-ray rejection, sky subtraction, and drizzling onto 
an output pixel scale of 0.05" . The narrowband and broad- 
band images of the same cluster were aligned using the lo- 
cation of ~ 20 bright stars. A narrowband excess image was 
constructed by direct pixel-to-pixel subtraction between the 
narrowband and broadband images, including an arbitrary 
scaling factor. We calibrated this scaling factor by checking 
that all bright cluster members, which are featureless ellip- 
tical galaxies with no emission lines in the respective filters, 
became consistent with the background in the excess image. 
For Abell 773 and Abell 68, the broadband images available 
in the archive did not directly overlap the Ha emission line, 
so an estimate of the broadband continuum was made by lin- 
ear interpolation between the adjacent F110W and F160W 
filters. 

The flux calibration of each image was verified using 
2MASS stars in the fields, and in all cases was found to 
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Figure 2. Ha excess images in the image plane (left) and reconstructed in the source plane (right). The image scales are in arcseconds 
in the image plane and in kpc in the source plane. Identified clumps are indicated in the source-plane images by black crosses, and the 
magenta ellipse shows the FWHM of the effective source-plane PSF, as described in the text. 



agree to within 15%, which is sufficient precision for our 
purposes. 

Colour HST images of the clusters are shown in Figure 
[1] with the critical lines at the redshift of the target arc over- 
laid. We use the transformation between image and source 
plane mapping from the best-fit cluster mass models (for 
details of t he mass models, see reference s in Table [1} with 
LENSTOOL l|Kneibll 19931 ; ljullo et aI1l2007t ) to reconstruct the 
images in the source plane, and show these in Figure [2] In 
order to reconstruct the source plane morphology, LENSTOOL 
uses the mapping between the image and source planes on 



a cluster-by-cluster basis and ray-traces the galaxy image. 
The lensing effect is to stretch the galaxy image - in most 
cases along one direction - and so the reconstruction cannot 
'create' new Hn regions, but rather the lensing has acted 
to extend them. As surface brightness is conserved by lens- 
ing, we then apply this conservation to obtain the intrinsic 
source plane flux. The total magnification is then simply the 
ratio of the image- to source-plane flux. To obtain the errors 
on the magnification, we use the family of best fit lens mod- 
els which adequately describe the cluster potential, derived 
by sampling the posterior probability distribution of each 
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parameter of the model (see I Richard et al.l (|201dl ) for more 
details). For each acceptable lens model, we reconstruct the 
arc and remeasure the amplification. We give the resulting 
magnification factors, /i, and associated errors in Table [1] 

In cases where the target is multiply-imaged, the im- 
ages were reconstructed separately and then adjusted for 
small differences in position and orientation before being 
combined. For MACS J0159, which consists of five images, 
only the first three were used due to the high magnification 
gradients in the fourth and fifth images resulting in high 
distortion in the source plane reconstructions. In the case 
of Abell 611, we use only the northernmost arc due to high 
distortion by a foreground galaxy lying close to the line of 
sight of the southern arc. 

We derive total magnification factors by comparing 
the total luminosities of the image- and source-plane Ha 
excess images. The intrinsic Ha luminosities are in the 
range 0.45 — 15 x 10 41 ergs~ 1 corresponding to SFRs of 
0.4 — 12MQyr _1 . These are at the faint end of the Ha lu- 
minosity function for this redshift range (see Fi gure EH), and 
probe fainter galaxies than the z ~ 2 sample of Ijones et al.l 
lj201Cfl . which covers the range 2.5 - 32 x 10 41 ergs"\ al- 
though the two samples overlap in luminosity. Due to the 
increased sensitivity provided by the lensing magnification, 
both of the lensed samples cover a lower range of intrinsic 
Ha luminosities than the sample of SINS galaxies studied 
by iForster Schreiber et al.l (|201lh , which were selected to 
have bright Ha and lie in the range 28 — 43 x 10 41 ergs _1 , 
making them rare, intensely star-forming galaxies. Thus, by 
harnessing gravitational lensing we are able to probe the 
more 'normal' star-forming population. 

Since gravitational lensing can preferentially shear one 
direction, we estimate the effective source-plane resolution 
by reconstructing the image of a star from the field reposi- 
tioned to lie at the centre of the target. The maximum lin- 
ear resolution, derived from the FWHM of the reconstructed 
star in the direction of greatest magnification, is 68-615pc 
with a median of 360pc, sufficient to resolve giant Hll re- 
gions. 

2.1 Comparison samples 

In order to interpret our high-z data, we exploit the Ha nar- 
rowband imaging f rom the Spitzer Infrar ed Nearby Galax- 
ies Survey (SINGS. iKennicutt et al]|2003h . which comprises 
Ha imaging of 75 galaxies with corrected SFRs of up 
to llMQyr -1 . We use the publicly available continuum- 
subtracted Ha narrowband imaging and restrict the sample 
to those with Ha detections with signal-to-noise of > 5 that 
have no significant defects in the galaxy images (determined 
by visual inspection). This restricts the SINGS sample to 41 
galaxies with SFR > 4 X 10" 4 M Q yr _1 . 

To ensure a fair comparison, we rebin the SINGS images 
so that the resolution is comparable to the high-z data and 
then threshold to the median surface brightness limit of the 
z ~ 1 — 1.5 observations. It is worth noting that thresholding 
the images in this manner excludes 10 — 50% of the total star 
formation. This should not affect the comparison between 
samples which have the same surface brightness- limit, but 
may serve as an indication of the fraction of star formation 
missed in high-z observations. 

To provide a comparison to local galaxies which are 
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Figure 3. Intrinsic Ha luminosities of the high-z s amples com- 
pared to Ha luminosity functions from HiZELS jSobral et al.l 
201lf). Also shown is the r ange of Ha luminosities of the 
Forster Schreiber et al. I fed 111) sample from the SINS survey at 
z ~ 2. The two lensed samples overlap in luminosity and are 
both at the faint end of the luminosity function, with the me- 
dian of the z ~ 1 — 1.5 WFC3 sample lower than that of the 
z = 1.6 — 2.6 OSIRIS sample by a factor of 6.6x, while the un- 
lensed SINS galaxies cover a range of higher Ha luminosities. 



more actively star-f orming, we use the VIMOS H a imag- 
ing spectroscopy of I Rodriguez- Z aurfn et all l|201ll) , which 
includes 38 LIRGs and ULIRGs at z < 0.13 with spatial 
resolution of 130 pc - 1.2 kpc and SFR < 25M yr" 1 . 

We also co mpare the z ~ 1 — 1.5 sample to the z ~ 
2 lensed arcs of Q ones et all lj2010T ). which were observed 
with Keck/OSIRIS. In order to provide a fair comparison, 
we have constructed narrowband images by summing the 
OSIRIS cubes over 100 A either side of the redshifted Ha 
emission line, matching the width of the WFC3 narrow- 
band filters. The resulting images are then corrected for 
lensing using the same image-to-source plane mapping as 
IJones et al.l ( 20101 ') in order to obtain the intrinsic galaxy 
properties. 

2.2 Determination of galaxy properties 

The total Ha luminosities of the galaxies in all samples are 
determined by summing all pixels in sky-subtracted images 
with signal-to-noise of > 3. In the case of the SINGS galax- 
ies, each image was checked by visual inspection and any 
foreground sources and defects masked. The resulting lu- 
minosities were then compared to the published values and 

found to agree to within ~ 20%. 

We convert Ha luminosity to S FR using t he Kennicuttl 
(|l998ah prescription, corrected to a IChabrierl (|2003T ) IMF. 



which reduces the SFR by a factor of 1.7 X. As we do 
not have constraints on the dust extinction, we adopt an 
estimate of Ana = 1 in all samples. This assumption is 
widely used in the literature a l thoug h it is the subject of 
some disagreement. [(Jam et all (|201Cf ) suggest a luminosity- 
dependent An a is more appropriate; were we to adopt 
their relation, we would obtain An a = 0.7 — 1.6 with a 
median Au a = 1.15. However , we also note that recent 
work bv lDomfnguez et al.l (|2012l ) suggests that galaxies with 
Lho < 4 x 10 41 ergs _1 may be consistent with having 
Au a = 0, and that above this threshold extinction increases 
in a luminosity-dependent way. Had we adopted this correc- 
tion instead, the SFRs of the majority of our galaxies would 
be reduced by a factor of 2.5 x . The exceptions are the three 
brightest z ~ 2 galaxies, in which the SFRs would increase 
by factors of 1.3 — 1.8x, and the z ~ 1 galaxies Abell 68 
and Abell 773; the former would be a factor of 1.8x lower, 
while the latter would be unchanged. Qualitatively, there 
is no significant impact on our results, as adopting either 
luminosity-dependent extinction relation would serve to in- 
crease the evolution we observe in H3.2.2I For simplicity and 
reproducability, we adopt Ana = 1 throughout. 

We define the sizes of the galaxies as twice the half-light 
radius. The half-light radius is determined using the contin- 
uum images to find the shape (i.e. the centre and major 
to minor axis ratio of an ellipse that best fits the galaxy), 
and then adjusting the semi-major axis of the ellipse until 
it encompasses half of the total Ha luminosity calculated 
in the manner described above. The galaxy- averaged star 
formation surface density, Esfr is defined from the total lu- 
minosity enclosed within two half-light radii per unit area. 



3 RESULTS AND ANALYSIS 

3.1 The Spatial Distribution of Star Formation 

A common theme in the recent literature is that high 
redshift galaxies are "clumpier" than galaxies in the 
local universe. This concept originates from the fre- 
quent appearance of "chain" galaxies in the high red- 
shift universe fe.g.lCowie et al.lll995l ; lElmegreen et al.ll2004l ; 



lElmegreen fc Elmegreenll2005h . Even without looking at the 
properties of individual star forming regions, it is interest- 
ing to compare the morphologies of the star-forming regions 
across the samples. 

From visual inspection, it is clear that there are sig- 
nificant differences between the samples. In particular, the 
surface brightness distributions of the galaxies show distinct 
differences in the different samples. In Fig. 2] we show the 
fraction of star formation in pixels above a given Esfr for 
the z ~ 1 — 1.5 and z ~ 2 samples, with the interquartile 
range of the thresholded SINGS sample shown for compari- 
son. 

To allow for the differing surface brightness limits of 
the samples, we only show star formation above a surface 
brightness of Esfr = 0.001 Mq yr _1 kpc -2 . This enables us 
to compare the star formation occurring in bright regions 
in a consistent manner. From the peaks - i.e. the points at 
which the lines tend to zero - we can see that the samples 
are different, with the z ~ 2 galaxies having peak surface 
brightnesses of around an order of magnitude higher than 
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the lower- z samples. Similarly, the z ~ 1 — 1.5 sample is sys- 
tematically brighter than SINGS sample, with the exception 
of MACS J1133, which is similar to the fainter z — galax- 
ies, MACS J0947 which is similar to the median of the z = 
sample, and Abell 773 which appears similar to the z ~ 2 
galaxies. 

As a statistical measure of the dumpiness of galaxies, 
we investigate using the Gini coefficient, G, which is used 
in econo mics to me asure the inequality of wealth in a pop- 
ulation It has values from to 1, where at 
the extremes G = for a completely uniform distribution, 
and G = 1 if there is only o ne non-zero value. Following 
iForster Schreiber et all (|201lh , we use it to quantify the dis- 
tribution of flux in an image, so a value close to one indicates 
that the profile has a single peak (in the case of G = 1, all 
of the flux would be in a single pixel), a galaxy with mul- 
tiple clumps would have a lower G, and at the extreme, 
a galaxy with completely uniform surface brightness would 
have G = 0. 

In the z ~ 1 — 1.5 sample, we find a narrow range of 
0.25 < G < 0.39 with a median of G = 0.34. The z ~ 2 sam- 
ple is marginally higher, with 0.42 ^ G ^ 0.56 and a median 
of G = 0.43. The 2 = SINGS sample has a similar median 
G — 0.45 but a much wider range of 0.05 G ^ 0.82, and 
the low-z ULIRGs have 0.38 < G < 0.85 with the high- 
est median G = 0.70. On the basis of the Gini coefficient 
there are no clear differences between the samples. Compar- 
ing the Gini coefficients with the visual appearance of the 
galaxies, the lack of distinction reflects the fact that a low 
Gini coefficient may arise from either a smooth distribution 
of star formation or from star formation that is concentrated 
into a large number of distinct clumps. Furthermore, we find 
no strong correlations between G and any of the properties 
of the galaxies. Clearly, to progress further we will have to 
compare the properties of individual clumps. In particular, 
we will show that the clump luminosity function provides a 
good means of distinguishing different galaxy star formation 
morphologies. 

3.2 Properties of star-forming clumps 

3.2,1 Definition of Clumps 

Studies of Hll regions or star-forming clumps have used 
a variety of methods to define and separate clumps 
from the background emission of the galaxy. Usually an 
isophote is defined at 3(7 above t he background no ise (e.g. 
iGonzalez Delgado fc Perez! Il997l ; IJones et al.1 1201(f ). How- 
ever, this method is clearly dependent on the noise proper- 
ties of the image, and thus is problematic when comparing 
local and high-redshift observations. In particular, as high- 
redshift galaxy images tend to have high relative noise levels 
and low dynamic range, the choice of isophote tends to se- 
lect only the brightest regions in the galaxy, neglecting any 
lower-surface brightness clumps and underestimating their 
sizes. 

An alternative is the IRAF task daof ind as employed 
bv lForster Schreiber et all (|201ll ). which is designed to lo- 
cate point sources in images. However, we found that it did 
not perform well on our sample. This is likely to be because 
daof ind requ ires an expected size of f e ature s to look for. As 
the clumps of lForster Schreiber et all |201ll ) are largely un- 



resolved, they were able to use the PSF of their observations 
as the expected size. As our clumps are resolved, the routine 
does not work reliably. Fo r this paper, we ther efore use the 
2D version of clumpf ind (| Williams et al.lll994T l. which uses 
multiple isophotes to define clumps. We defined the contour 
levels with respect to the rms noise in the image, starting 
at 3a and increasing in la intervals until the peak value 
of the image is reached. The data are first contoured at the 
highest level to locate clumps, and the algorithm then works 
down in brightness through the contour levels. Any isolated 
contours are defined as new clumps, while others extend ex- 
isting clumps. If a contour surrounds one existing peak, they 
are allocated to that clump, and any which enclose two or 
more are divided using a 'friends-of-friends' algorithm. The 
advantages of this approach are that it enables a consistent 
clump definition to be applied to multiple data sets, lower 
surface brightness clumps are not excluded, and there is no 
assumption made about the clump profile. 

The clumps identified by clumpf ind were all confirmed 
by visual inspection to remove any sources not associated 
with the target galaxy, of particular importance in the case 
of the SINGS images where foreground sources lie close to or 
overlap the target galaxies. The area A of the clump is then 
obtained from the number of pixels assigned to it, multiplied 
by the source-plane pixel scale, and from this we define the 
effective radius r — yjA/it. We only accept clumps where 
2r is larger than the FWHM of the PSF, so all clumps are 
resolved. 

Due to the manner in which clumps are 'grown,' their 
sizes returned by clumpf ind tend to be larger than those 
obtained by other methods. As a comparison, we also fit 
a 2D elliptical Gaussian profile to each peak and mea- 
sure the FWHM. A comparison of the clump radii found 
by the two methods is shown in Figure [5] The rms dif- 
ference between the two radii is ~ lOOpc, and on aver- 
age we find that clumpf ind o utput s sizes 25% higher than 
the FWHM. IWisnioski et all (|2012T l note that sizes defined 
through isophotes can be unreliable due to the level of 'tun- 
ing' required to select an appropriate isophote level in a 
given galaxy. This is less significant with clumpf ind because 
this tuning is not required; the use of multiple isophote lev- 
els in all galaxies allows the levels to be defined in a consis- 
tent way across a large sample. We therefore find much lower 
scatter between the isophotal sizes output by clumpf ind and 
the clump FWHM than they do in their sample. Through- 
out this work, we use the clumpf ind size for all samples, and 
give error bars that encompass the FWHM of the clumps. 

The sizes and Ha-derived SFRs of the z ~ 1 — 1.5 clumps 
are given in Table [2] We analyse these properties in compar- 
ison to the other samples below. 

3.2.2 Clump properties 

One way of quantifying the 'dumpiness' of a galaxy is to 
consider the fraction of a galaxy's total Ha luminosity con- 
tained within clumps. We find medians of 31% in SINGS, 
36% for the z < 0.13 ULIRGs, 50% for the z ~ 1 - 1.5 sam- 
ple and 68% for the z ~ 2 sample. Thus, as expected, the 
higher- z galaxies are clumpier than their local counterparts. 

We now consider the properties of the clumps them- 
selves, and first compare the Ha-derived SFR to the clump 
radius, as shown in Figure [6] Locally, there is a well-defined 
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Figure 4. The fraction of star formation within each galaxy occurring above a given surface brightness, for the z ~ 1 — 1.5 and z ~ 2 
samples. The shaded region is the interquartile range of the SINGS z ~ sample. There are two galaxies, MACS J1133 and MACS J0947, 
from the z ~ 1 — 1.5 sample with similar surface brightnesses to the z = sample, and the remainder are systematically brighter. The 
z ~ 2 sample has significantly higher surface brightnesses. Hence, there is clear evolution in the surface brightnesses of galaxies with 
redshift. 



relatio nship between these properties, as found bv lKennicuttl 
(| 19881 ) who found almost constant surface brightness in 
local Hll regions, except in merging and interacting sys 



tems l|Bastian et al] 2006|). The situation at hig h-z, though, 



appea rs different; Swinbank et al.l (|2009l ) and Ijones et al.l 
lj2010T l found clumps with SFRs of ~ 100 x higher at a given 
size than found locally, in systems with no evidence of in- 
teractions. 

Figure [ 6] is a n updated version of one presented in 
Ijones et all (|201fJ ). where we have re-analysed the z ~ 2 
and SINGS galaxies using clumpf ind so that clumps are de- 
fined consistently across all samples, and we have added the 
results from our new z ~ 1 — 1.5 data set and the z < 0.13 
ULIRGS as well as the z = 1 - 2 results from SHlZELS 
(Swin bank et al. in prep.) and WiggleZ (|Wisnioski et al.l 
2012). We show lines of median surface brightness in the 



samples, and vertical offsets from these lines represent dif- 
ferences in the surface density of star formation, Esfr, in 
the clumps. We will explore the relation of these offsets to 
global galaxy properties in Section [4] 

We note that the clumps we identify in the SINGS 
galaxies are derived from images which have been degraded 
to comparable resolution to the high-z data, and we find the 
effect of this is to decrease the surface brightness by a fac- 
tor of ~ 2 x , as the size increases more than the luminosity. 
The points in Figure [6] move along the vector labelled 'A'. 
Defining clumps in the z — sample in this way ensures the 
fairest possible comparison with the high-z data. 

Upon re-analysis using clumpf ind, we find some lower 
Ssfr regions in the Ijones et al.l (|2010l ) sample, but they 
all remain separated from the local relation by a factor 
of ~ 100 x. This confirms the large differences between 
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Table 2. Properties of clumps identified in the z < 
determined as described in £13-2,11 



1 — 1.5 sample, 



Clump 



radius (pc) SFR (M yr- 1 ) 



MACS J0947-1 


350 ± 56 


0.054 


± 





010 


MACS J0947-2 


324 ± 22 


0.0328 


± 





0086 


MACS J0947-3 


384 ± 48 


0.045 


± 





012 


MACS J0947-4 


334 ± 39 


0.0350 


± 





0092 


MACS J0947-5 


318 ± 38 


0.0339 


± 





0089 


MACS J0947-6 


376 ± 6 


0.0347 


± 





0091 


MACS J0947-7 


311 ± 17 


0.0261 


± 





0068 


MACS J0947-8 


149 ± 9 


0.0050 


± 





0013 


MACS J0159-1 


402 ± 89 


0.282 


± 





060 


MACS J0159-2 


370 ± 77 


0.203 


± 





043 


MACS J0159-3 


530 ± 130 


0.355 


± 





076 


MACS J0159-4 


468 ± 13 


0.170 


± 





036 


Abcll 611-1 


730 ± 180 


0.370 


± 





083 


Abell 611-2 


560 ± 160 


0.181 


± 





041 


Abcll 611-3 


630 ± 140 


0.200 


± 





045 


Abell 611-4 


390 ± 59 


0.065 


± 





015 


Abell 68-1 


378 ± 26 


0.081 


± 





016 


Abell 68-2 


132 ± 8 


0.0075 


± 





0015 


Abell 68-3 


375 ± 24 


0.076 


± 





015 


Abell 68-4 


354 ± 25 


0.070 


± 





014 


Abell 68-5 


509 ± 31 


0.114 


± 





023 


Abell 68-6 


337 ± 33 


0.062 


± 





013 


Abell 68-7 


386 ± 61 


0.099 


± 





020 


Abcll 68-8 


205 ± 42 


0.0273 


± 





0055 


Abell 68-9 


348 ± 44 


0.069 


± 





014 


Abcll 68-10 


299 ± 89 


0.066 


± 





013 


Abell 68-11 


328 ± 31 


0.057 


± 





012 


Abell 68-12 


312 ± 15 


0.0476 


± 





0097 


Abell 68-13 


412 ± 60 


0.100 


± 





020 


Abell 68-14 


293 ± 16 


0.0429 


± 





0087 


Abell 68-15 


263 ± 5 


0.0302 


± 





0061 


Abcll 68-16 


280 ± 50 


0.0454 


± 





0092 


Abell 68-17 


169 ± 7 


0.0132 


± 





0027 


Abcll 68-18 


195 ± 33 


0.0215 


± 





0044 


Abell 68-19 


224 ± 6 


0.0199 


± 





0040 


Abcll 68-20 


239 ± 51 


0.0344 


± 





0070 


Abcll 68-21 


135 ± 5 


0.0069 


± 





0014 


Abell 68-22 


163 ± 23 


0.0142 


± 





0029 


Abcll 68-23 


112 ± 2 


0.00489 


± 





00099 


Abell 68-24 


171 ± 31 


0.0161 


± 





0033 


Abell 68-25 


98 ± 5 


0.00416 


± 





00085 


Abcll 68-26 


122 ± 8 


0.0066 


± 





0013 


Abcll 2390-1 


352 ± 59 


0.115 


± 





025 


Abcll 2390-2 


404 ± 68 


0.136 


± 





030 


Abcll 2390-3 


409 ± 6 


0.086 


± 





019 


Abcll 2390-4 


366 ± 26 


0.067 


± 





015 


Abcll 2390-5 


463 ± 9 


0.111 


± 





024 


Abcll 2390-6 


347 ± 13 


0.059 


± 





013 


Abell 2390-7 


470 ± 1 


0.093 


± 





020 


Abcll 2390-8 


341 ± 73 


0.069 


± 





015 


Abcll 773-1 


1040 ± 200 


5.6 


± 


1 


3 


Abcll 773-2 


1430 ± 180 


9.6 


± 


2 


2 


MACS J1133-1 


1120 ± 100 


0.118 


± 





025 


MACS J1133-2 


890 ± 160 


0.068 


± 





015 


MACS J1133-3 


790 ± 280 


0.057 


± 





012 


MACS J1133-4 


835 ± 4 


0.0334 


± 





0072 


MACS J1149-1 


174 ± 34 


0.084 


± 





020 



1000 



o 
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Figure 5. Comparison of the clump size r c i ump fi n( j output by 
clumpfind with the size t"fwhm obtained by taking the FWHM of 
a 2D Gaussian profile fit. On average, clumpfind outputs sizes 25% 
larger than the FWHM. For consistency, we adopt the isophotal 
size output by r c i um pfi n d in all samples. 



the local and high re dshi ft population alrea dy noted by 
ISwinbank et all (|2009l) and I Jones et alj (^Olfl ). 

Our new z ~ 1 — 1.5 sample fits in between the SINGS 
and z ~ 2 samples, with the exception of the two regions 
from the most compact source Abell 773, which have Esfr 
similar to the z ~ 2 sample, and the four regions from 
MACS 1133, which are similar to z = clumps. This in- 
dicates clear evolution in clump surface brightness, Espr, 
with redshift. 

The surface brightness limit of the z ~ 2 data means 
that we cannot identify the low star formation rate clumps 
in that sample. We show a dotted line representing the lower 
limit at which we define clumps in the z ~ 2 galaxies. It is 
likely that there are additional clumps which he below this 
limit and are undetected; however, such clumps make only 
a small contribution to the total SFR, as we shall discuss in 
Section 13.31 

Selection effects have no impact on the lack of high sur- 
face brightness regions in the lower redshift samples, how- 
ever. The intense star-forming regions are clearly more com- 
mon in high- 2 galaxies; they are found only in extreme sys- 
tems such as the Antennae locally, but exist in all five of the 
2 ~ 2 galaxies and one of the eight z ~ 1 — 1.5 sample. 

As noted in § [2] the z ~ 2 galaxies have 'normal' SFRs 
for their redshift, below the knee of the Ha luminosity func- 
tion. The offsets seen in the figure emphasise the impor- 
tance of analysing clumps in terms of their surface bright- 
ness. This is even more evident if the clumps belonging to a 
single galaxies are examined separately. Rather than being 
distributed across the plot at random, individual galaxies 
form a much tighter sequence with all the clumps sharing a 
common surface brightness, particularly in the low redshift 
sample. Thus the spread in clump properties in Figure [6] 
appears to be driven by global differences in the galaxies. 

We therefore next compare the clump Espr to the prop- 
erties of their host galaxies in Figure [7] In the left-hand 
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Figure 6. Ho SFR for extracted Hn regions as a function of size, co mpared to the lensed z ~ 2 sample of I Jones et al.l ll2010h. high-z un- 
lensed samples from SHiZELS (Swinba nk et al. in prep) and W iggleZ l lWisnioski et al.l2012t) , low-z ULIRGs from Rodrfguez -Zaurfn et al,l 
ll201ll ) and the z = SINGS galaxies l lKennicutt et alj|2003h . Star formation rates are calculated using the iKennicuttl dl998afT prescrip- 
tion adjusted for a Chabrier IMF with a dust extinction Ajj a = 1 in all samples, and the error bars of the high-z lensed sources are 
dominated by the uncertainty in the lensing magnification. Dashed lines show the median surface brightnesses in the SINGS, z ~ 1 — 1.5 
and z ~ 2 samples. The black dotted line indicates the sensitivity limit of the z ~ 2 OSIRIS observations. The arrow indicates the effect 
of degrading the image resolution, as discussed in the text. The four lowest surface brightness clumps in the z ~ 1 — 1.5 sample come from 
one galaxy (MACS J1133), and the two brightest regions arc from Abcll 773, the most compact galaxy in the sample. The remaining 
galaxies have clumps with surface brightnesses in between those of the z = and z ~ 2 samples, similar to local ULIRGs. 



panel, we correlate the clump properties with the total star 
formation rate of the galaxy. For clarity, we plot the median 
clump Ssfr in each individual galaxy, and the error bars 
encompass the central 68% of clumps within each galaxy 
(i.e. 1 a if they follow a Gaussian distribution) There is 
evidence for correlation between the clump Esfr and the 
galaxy Ha luminosity (which we assume to be proportional 
to the total SFR); we find a Spearman rank correlation coef- 
ficient p = 0.69, representing a 5.8 a deviation from the null 
hypothesis of no correlation. This suggests that the star- 
formation in the high-z sample follows a similar trend to 
the local sample, and that the differences seen in Fig. [6] may 
arise from the higher total star formation rates of the high 
redshift galaxies. 

For the majority of the samples, an even stronger re- 
lation arises if we compare the clump Esfr to the galaxy- 
averaged Ssfr. This is shown in the right-hand panel of 
Figure and has a correlation coeficient p = 0.79 with 
6.6 a significance. The ratio of clump-to-average Esfr can 
be thought of as a measure of the 'dumpiness' of the galaxy. 



We conclude that the properties of star-forming clumps 
in a galaxy are strongly dependent on the global Esfr of the 
galaxy. Galaxies with higher overall Esfr have higher clump 
surface densities and are correpondingly offset in the clump 
size - star formation rate relation. While this accounts for 
some of the differences seen in Fig. [6l it is also clear that 
there are more bright clumps in the higher redshift galaxies. 
We quantify this below. 



3.3 Hn region luminosity functions 

A quantitative measure of the clump brightness is to con- 
struct a luminosity function (LF) of Hh regions. In the loca l 
univ erse, the Hn LF is presente d in Kennicutt et al.l (|l989l ) 
and iGonzalez Delgado fc Perezl (1 19971 ). They demonstrate 
that the LF can be fitted by a broken power-law, or by a 
power-law with an exponential break. In order to be con- 
sistent with our definitions of clump sizes, we reanalyse the 
local data in order to construct our own LF. The results 
are shown in Figure [8] The left panel of the figure shows 



© 2011 RAS, MNRAS 000, [THS] 



Evolution in the size and luminosity of giant Hll regions 11 




0.0010 - 



0.0001 I I I I I I I I_ 

10 39 1 40 1 4 ' 10" 2 0.01 0.10 1.00 

Ha Luminosity (erg/s) Galaxy average S SFR (M yr~'kpc~ 2 ) 



Figure 7. Comparisons between the star formation surface density EgpR of star-forming clumps within each galaxy and the intrinsic 
Ha luminosity and galaxy-averaged EgpR. The clump SgpR shown is the median for each galaxy, with error bars encompassing the full 
range of Ssfr f° r an clumps within each galaxy. The solid line is the best fit to the data, and the dashed line illustrates the clump Esfr 
expected from theory, discussed in Section [4] We find that both are correlated at the 5cr level, implying that we find more high-EgpR at 
high redshift because there are more high-SFR and Sgppi galaxies at this epoch. 
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Figure 8. Cumulative luminosity functions of Hll regions in the SINGS, z ~ 1 — 1.5 and z ~ 2 samples, shown as a mean per galaxy 
(left) and normalised by total galaxy Ha luminosity (right). We plot the SINGS sample as a whole and subd ivded into three luminosity 
bins. The shaded regions illustrate model predictions from the GMC mass functions of lHopkins et al~ | J2012T) for Milky Way-like (grey) 
and high- 2 (blue) simulations. We find evolution in the Hll region luminosity function with redshift, which seems to be driven by the gas 
fraction of the disk. 



the cumulative number of regions per galaxy as a function 
of Ha luminosity. The normalisation of each bin takes into 
account the different surface brightness limit of the galaxies, 
with error bars computed from the Poisson error in counting 
regions. The slope of the power-law part of the mass function 
is ~ —0.75, so that although the LF appears steep in this 



representation, most of the total luminosity is contributed 
by the brightest Hn regions. 

Solid black points show the average of all SINGS galax- 
ies. However, since we will be comparing the galaxies cover- 
ing a range of luminosities and redshifts, we have separated 
the galaxies from the SINGS sample into 3 total Ha lumi- 
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nosity bins. At a fixed luminosity, galaxies with lower total 
emission have fewer regions, but the shape of the luminosity 
function is similar. In order to emphasise the similarity of 
the mass function, we normalise each of the curves by the 
total star formation rate of the host galaxies, and show this 
in the right-hand panel. The similarity of all the Hn region 
LFs is now clear. 

There is a striking similarity between the LF of the 
z ~ 1 — 1.5 sample and that of the highest SFR galaxies in 
the low- 2 SINGS sample. The excess of very bright regions 
(L ~ 10 41 ergs -1 ) is down to one galaxy, Abell 773, which 
is the same compact galaxy for which we found the clump 
surface brightnesses to be more typical of the highest red- 
shift galaxies. The low luminosity slope of the LF tends also 
to be flatter than that seen in the low redshift galaxies, but 
it is hard to quantify this difference without directly com- 
parable surface brightness limits and is likely to be affected 
by unresolved regions which are excluded. In any case, these 
faint regions contribute little to the total flux. 

In both panels, the Hn region LF for the highest redshift 
galaxies is strongly offset from the relation seen in the low 
redshift SINGS sample and from the sample at z ~ 1 — 
1.5, but is similar to the low-redshift ULIRGs. Although the 
data do not probe the low- luminosity slope of the LF, these 
galaxies have much brighter regions than are seen at lower 
redshift. The right-hand panel emphasises that the this is 
not because they contain many more regions overall. 

In order to compare our data to models of mass func- 
tions, we must relate the measured Ha luminosities to model 
clump mass, M. As an estimate, we use the Ha-derived 
SFR and adopt SFR (A^yr -1 ) = 4.6 ± 2.6 x 1O _8 M 
|Lada et al.ll201oT ). This empirical relation is based on local 
molecular clouds and applies to the high-density gas where 
Ak > 0.8 mag. However, we note that this relation is con- 
sistent with the far-infrared-derived star formation rate and 
CO-derived gas masses of star-forming clumps r eported in a 
lensed z = 2.3 galaxy bv lSwinbank et all (|201ll 'l. but clearly 
more high-resolution CO observations of high- 2 galaxies are 
required to confirm this. As a guide, we include this conver- 
sion on the upper axis of Figure [5] 

The shapes of the LFs can be approximated by a power 
law with an exponential cut-off at some high luminosity or 
mass. The difference between the samples' LFs is then best 
described by a pure luminosity evolution, so that the cut-off 
shifts to higher luminosity/mass at higher redshift. 

To demonstrate this, we include shaded regions repre- 
senting a Schechter function of the form 

where we adopt the median value of a = —0.75 from 
iHopkins et"all <|2012h . The normalisation iVo is arbitrary, so 
we fit No to the z = data and then keep it fixed while 
allowing Mo to vary in order to find best-fit values of in the 
different samples. The best-fit values are Mo = 4.6I2J x 
10 7 Mq at z = 0, M = S.Ot 1 ^ X 10 7 Mq at z ~ 1-1.5 and 
Mo = l.St 10 x 10 9 Mq at z ~ 2, where the errors are found 
with a bootstrap method. We shade the best-fit Schechter 
functions at z = and z ~ 2 in grey and blue respectively 
in Figure [8] The normalisation of the model in the right- 
hand panel is obtaines by summing the luminosities of the 



clumps. The result is in remarkably good agreement with 
our observations. 

Not only do the highest-redshift galaxies have Hn re- 
gions that are higher surface brightness, but the character- 
istic luminosity of the regions is higher too. We suggest that 
the presence of high-luminosity regions may be a good op- 
erational definition of the dumpiness of a galaxy. 



4 DISCUSSION 

In the previous section, we presented an analysis of star 
forming regions in galaxies at 2 = 0, 2 ~ 1 — 1.5 and z ~ 2. 
We find that the luminosities of the regions in z = 1 galaxies 
are similar to those of bright (L > 10 40 ergs _1 ) galaxies at 
low redshift, but the surface brightnesses are systematically 
higher. At higher redshifts, the properties of the galaxies 
change, with the galaxies having clumps that are both much 
higher surface brightness and shifted to higher total lumi- 
nosities. This accounts for the qualitative impression that 
the most distant galaxies are "clumpier." 

We also noted that the increase in the surface brightness 
of Hn regions tracks the increase in the average star forma- 
tion rate surface density, Y*sfr, of the galaxies. The obser- 
vations are consistent with the changing properties of the 
Hn regions being driven by changes in the overall T,sfr of 
the galaxies. We can link the increase in the observed T,sfr 
to an increase in the gas surface density of these galaxies 
by assuming that the Kennicutt-Schmidt law holds at 2 ~ 2 
as well as at z = 0. In this case, we have an emerging pic- 
ture that the changes we see are likely driven by greater gas 
surface densities at higher redshift. 

The connection between the increasing surface density 
of clumps and th e greater peak brightness arises naturally 
from this picture (|Hopkinsll201ll ). The clump mass required 
for collapse on scale 7? from a turbulent ISM is given by the 
Jean's mass, Mj: 

^=^ M ' M B^G^' (2) 
where at{R) is the line of sight turbulent velocity dispersion. 
Assuming a turbulent velocity power spectrum E(k), the 
velocity dispersion ot{R) for wavenumber k = 1/R is 

a\ = kE(k) oc R?' 1 , (3) 

wher e p m 2 for supersonic turbulence appropriate to the 
ISM (|Burgerdll97l ; ISchmidt et alj|2009h . 

To determine the normalisation of the relation, we as- 
sume that the clumps are located in a marginally unstable 
disk. We note that the available kinematic data for the 2 ~ 2 
sample and for MACS J 1149 supp ort this assumption, as do 
larger surveys (|Genzel et alj|201lh : nonetheless, clearly this 
is uncertain without dynamical data for the entire sample. 
However, if we make the assumption that the galaxies are 
rotating disks with Toomre parameter Q fa 1, we can relate 
the epicyclic frequency, n of the disk to its scale height, h: 
k w a t (h)/h. 

7rGEo 7rGEo/l 

Since the stability of the disk is a global phenomenon, we 
will associate Eo with the average surface density of the 
star forming disk, E c ii s k, and treat the quantities appearing 
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in Eq. [4] as appropriate global disk averages. Since the disk 
is made up of both stars and gas, we must take an appropri- 
ate avera ge of the surfac e densities in the two components. 
Following Rafikovl (|200l| ), and focussing on the largest un- 
stable fluctuations, the appropriate combination of gas and 
star surface densities (denoted E s and E.) is 



Edisk — Eg -f- 



1 + 



(5) 



where f a — cr* /at is the ratio of the velocity dispersion of 
the stellar component to that of the gas. 
Assuming Q ~ 1, we can write, 

a 2 t {R) = (7rGE disk /i) fjY « nGT, disk R, (6) 

where we have used Eq.[3]to relate a t (h) to a t (R). Combin- 
ing with Eq. [2] gives a critical density for collapse of 



Pc(R) 



8 disk iT 



Assuming that the cloud contracts by a factor 
collapses, the post-collapse surface density is 

Ecloud ^ 10p c R ~ 10 Edisk- 



(<) 

2.5 as it 

(8) 



Thus, for the turbulent power spectrum p w 2, the sur- 
face density of collapsed clouds is independent of radius and 
proportional to the surface density of the disk. The nor- 
malisation of the relation follows from the collapse factor 
and the requirement that the disk is marginally stable. This 
model provides a good description o f clouds in the Milky 
Way (Larson's laws) as discussed in IHopkinsl l |20Tlf ). and 
predicts that the surface brightnesses of clouds should in- 
crease as the gas surface density (and thus overall average 
star formation rate surface density) increases. If we assume 
a constant Kennic utt-Schmidt law of the form E gas oc Egp R 
l|Kennicuttlll998tJ ). we can compare Eq. [2] to our data; we 
therefore overplot this line on Figure [7] and find that it is in 
good agreement with the observations. 

Moreover, the model predicts that the most massive 
clouds should increase in size with the average gas surface 
density. This follows from the marginal stability condition 
(Eq. [4j , since density structures on scales greater than h will 
tend to be stabilised by disk rotation. 

This can be demonstrated formally by examin- 
ing the dispersion relation for a fin i te th ickness disk 
|Begelman fc Shlosmar] 120091 '). IHopkinsl (|201lh shows that 
this leads to an exponential cut-off of the clump mass func- 
tion above a mass 



M 



4-7T c 



(9) 



where we have used Eqs. [3] and [4] to express h as a 
function of Edisk and k = fdisk/-Rdisk (where Wdisk is the 
disk circular velocity and i?disk is half-mass radius of the 
disk). Expressing k in units of 100 km s" 1 kpc -1 we obtain 
a normalisation of 



Ml = s 6 x 10 3 I — 

M e ' X VlOM pc-V VlOOkm/skpc- 1 



(10) 

We can check that this results in a reasonable value 



of Mo in the Milky Way by using a gas surface den- 
sity E gas ~ 10 M(7) pc~ 2 and a g as fraction of 10% 
with fa- — 2 (Korchagi n et al] 120031 ) to obtain an effec- 
tive Edisk ~ 35M0pc~ 2 . With k = 36.7 km s" 1 kpc -1 
(|Feast fc Whitelocklll997h . this gives M ~ lO 7 Af , in good 
agreement with the character istic mass of the largest Galac- 
tic GMCs (|Stark fc Leell2006h . 

Eq. [10] shows that the mass cut-off depends strongly on 
the disk surface density — the higher the surface density, the 
more massive the clumps that are able to form. This trend 
can, however, be opposed by the stablising effect of angular 
rotational speed. For a fixed disk radius, a higher circular 
velocity tends to reduce the mass of the largest clumps. 

For low-redshift galaxies, simple theoretical models sug- 
gest that Rdisk oc ^disk since halo spin is weakly dependent 
on the halo mass (|Mo et alj|l998l ), and thus we should ex- 
pect the dependence on the disk surface density to be the 
dominant trend controlling the cut-off clump mass. This is 
confirmed by analysis of the obser ved properties of galaxies. 
For example, IDutton et ail (|201lT) find 



-Rdisk 

kpc 



2.5 



«disk 

100 km/s 



(11) 



in the local universe. Combining this with the observed de- 
pendence of the disk rotation velocity on galaxy mass 



10 10 A/ ( 



0.25 







/ «disk 

V100 km/s 



we obtain 



100 km/s kpc- 1 



0.38 



1O 1O M 



(12) 



(13) 



which shows that k is very weakly dependent on the galaxy 
mass, and the variation in the clump mass functions of local 
galaxies is driven is driven by the disk surface density. 

We also note that in a disk of constant circular veloc- 
ity, k scales with radial distance r as n oc r _1 , while the 



gas surface den sity profiles are shallow, with E gas oc 



-4/3 

(|Fu et alJl2O10l ). Thus from Eq. 1101 there is no dependence 
of Mo on r; while the surface density is higher towards the 
centre of the galaxy, this is balanced by the higher rotational 
frequency. This explains the observations that clump prop- 
erties appear to be driven by the global properties of their 
host galaxies rather than by local conditions, and this allows 
us to use disk-averaged values of k and E. 

We have no measurements of the gas contents of our 
samples, but dynamical information available for the z ~ 2 
sample permits us to predict the cut-off mass in these galax- 
ies from our model if we estimate Edisk from the dynami- 
cal m ass. Using the measurements reported in Ijones et al.l 
l|2010l . see their Table 2), we compute a cut-off mass ranging 
from 3.3 x 10 6 - 3.1 x 10 9 M for the z ~ 2 sample. The me- 
dian value is 5 x 10 7 Mq , approximately 5 x higher than the 
Milky Way value. We therefore expect that the z ~ 2 sam- 
ple should contain clumps of higher mass and luminosity, as 
observed. However, the uncertainty in the cut-off mass for 
the z ~ 2 sample is very large due to uncertainties in Edisk 
and ft, which prevents us from making a precise comparison 
of the cut-off mass in the different samples. 

To understand how the dumpiness of galaxies evolves, 
we must therefore use simula tions to estimate th e evolu- 
tion of their scaling relations. IDutton et al.l (|201ll ) present 
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Figure 9. The evolution of cut-off mass Mq and clump surface density S c i ump in comparison with model predictions (equation 1 10| l. The 
model is strongly dependent on the assumed evolution of gas fraction with redshift; we show tracks for / gas oc (1 + 2) 2 ±0.5 anc j ^ ne t n j c k 
arrows show the effect of increasing the gas fraction. We also show a track for / CT = 1 - i.e. assuming that the gas and stars have the same 
velocity dispersion. This makes the disk unstable on larger scales and would lead to higher-mass clumps at low-z than are observed. The 
impact of this is reduced at higher-z where gas dominates the disk dynamics. As increases, the points move in the direction of the 
arrow in the lower right corner. The model provides a good fit to the data, demonstrating that larger, higher surface brightness clumps 
at high-z are a natural consequence of increasing gas fractions, which explains the observed morphologies of high-z galaxies. 



a simple analyti c model that seems to describe the observa- 
tional data well dTruiillo et al.ll2006l ; iForster Schreiber et alj 



|2009| ; I Williams et al-lboipfT We use their scaling relations for 
mass, size and rotational frequency with redshift in com- 
bination with Eq. [S] and Eq. [TO] to predict how the cut- 
off mass and clump surface brightness should evolve with 
redshift. Figure [5] illustrates this evolu tion for a gas frac 



tion evolution of / gas oc (1 + z) 



(|Geach et al.ll201lD . 



The arrows show how altering the assumed gas fraction 
changes the model. This suggests that the changing clump 
properties are a natural consequence of increasing gas frac- 



tions dominating high-z galaxy dynamics. The high gas 
fractions probably arise from high gas infall rates at high 
redshift dBournaud fc Elmegreenl |2009| ; iKrumholz fc Dekel 
l20ld ; iBournaud et al.ll201lh ; however, our observations do 
not directly rule in or out cold flows. Our results merely 
require high gas fractions, and cold flows are a method of 
maintaining the gas supply. Crucially, we note that this ef- 
fect is tempered by the more compact nature of galaxies, 
which leads to higher epicyclic frequencies that limit the 
collapse on larger scales. The need to include the k term is 
apparent from our Hn region luminosity functions: without 
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it, a factor 10 increase in disk surface density would corre- 
spond to an increase in clump luminosity of 1000 x, and we 
do not observe such a large increase. 

To summarise, we find that our simple theoretical model 
is in good agreement with the observations and suggests that 
the evolving 'dumpiness' of galaxies is a manifestation of the 
different cut-off mass of the Hn region luminosity function, 
which is driven by evolution in the gas fraction with redshift. 



5 CONCLUSIONS 

We have used HST/WFC3 to obtain narrowband Ha imag- 
ing of eight gravitationally lensed galaxies at z ~ 1 — 1.5. 
The magnification provided by the lensing enables us to 
reach spatial resolutions in the source plane of 68 — 615 pc. 
In addition, to provide comparisons we have re-analysed 
the lensed z ~ 2 sam ple observed with Keck / OSIR IS by 
iJones et al.l <|2010h . the kodriguez-Zaurm et ail (|201ll ) sam- 
ple of z < 0.13 (U)LIRGs observed with VLT/VIMOS and 
the Ha narrowband im aging of the z = SINGS survey 
l|Kennicutt et alj|2003h . 

The high-z samples have 'clumpy' morphologies, domi- 
nated by a few large regions of high Ha luminosity, which we 
use as a proxy for the SFR. We have extracted star-forming 
clumps from the galaxies in each sample and examined their 
properties. The clumps follow similar SFR-size scaling rela- 
tions in all samples, but the normalisation of the relation 
exhibits systematic offsets to higher surface brightness at 
higher redshifts. The normalisation appears to be approx- 
imately constant within a given galaxy, implying that this 
relation is driven by global galaxy properties. 

On comparison with the properties of the host galax- 
ies, we find that all samples follow approximately the same 
scaling relations between the clump surface brightness and 
both the host galaxy's total Ha luminosity, Lua , and its av- 
erage surface density of star formation, Esfr, and that they 
evolve along this relation in decreasing Lua and Esfr with 
decreasing redshift. 

We have measured the luminosity function of clumps in 
the samples, and shown that the z ~ 1 — 1.5 sample is similar 
to the higher-Z/Ha members of the SINGS sample. When 
normalised by the host galaxies' total SFR, the SINGS and 
z ~ 1 — 1.5 samples can be fit by the same Schechter function, 
while the ULIRGs and z ~ 2 samples are offset horizontally. 
This shift can be explained by an increase in the cut-off 
mass of the Hll region luminosity functions of the ULIRGs 
and z ~ 2 disks. 

We present a simple theoretical model which shows that 
the evolution in luminosity and surface brightness are con- 
nected, and are driven by the competing effects of disk sur- 
face density Edisk an d the epicyclic frequency re. Galaxies 
at high redshift tend to have higher Edisk, which increases 
the maximum mass of clumps that are able to form; how- 
ever, this is tempered by the more compact nature of high-z 
galaxies, implying higher re, which impedes collapse on the 
largest scales. 

We have shown that this model is consistent with the 
evolution in clump properties seen in our data. We there- 
fore conclude that the clumps observed in high-z galaxies 
are star-forming regions analogous to those found locally 
but with higher masses and surface brightnesses. As Hll re- 



gions in the distant Universe are larger and brighter, they 
give rise to the 'clumpy' appearance. The increase in clump 
luminosity is driven primarily by increasing gas fractions 
at high-z. This clearly motivates further study with ALMA 
to better quantify the evolution of gas properties in high-z 
galaxies. 
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